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ABSTRACT 

The Laplace operator admits infinite self-adjoint extensions when considered on a 
segment of the real line. They have different domains of essential self-adjointness char- 
acterized by a suitable set of boundary conditions on the wave functions. In this paper 
we show how to recover these extensions by studying the continuum limit of certain dis- 
cretized versions of the Laplace operator on a lattice. Associated to this limiting proce- 
dure, there is a renormalization flow in the finite dimensional parameter space describing 
the dicretized operators. This flow is shown to have infinite fixed points, corresponding 
to the self-adjoint extensions characterized by scale invariant boundary conditions. The 
other extensions are recovered by looking at the other trajectories of the flow. 



1. Introduction 



Since the time of their introduction by Wilson lattice field theories have proven to 
be a successfull tool in the study of quantum field theories. In this approach, one replaces 
the space-time continuum (or rather, in the Hamiltonian version of the method 0, the 
space continuum ) with a discrete lattice. Correspondingly, the action (or the Hamilto- 
nian) describing the continuum theory is replaced by a discretized one. The spacing a 
of the lattice then provides a natural cut-off for the theory and the continuum limit is 
recovered by taking a to zero. This process is governed by the renormalization group, 
which dictates how the parameters on which the discretized action depends have to evolve 
with a so that, for a sufficiently small, the physics becomes insensitive to the presence of 
the lattice spacing; when this happens, one says that the continuum limit exists. Unfor- 
tunatly, for most of the field theories of interest, the analysis of the renormalization group 
is too complicated to be carried out exactly and can be attacked only within perturbation 
theory. 

In this paper we discuss a non trivial illustration of the ideas of Wilson already in the 
context of quantum mechanics. The idea of renormalization has been applied already to 
quantum mechanics by other authors |Q, but in a different way from that pursued 
here. In particular, in 0], it is shown how the use of renormalization allows one to make 
physical sense of the inverse square potential, in one space dimension, by rendering finite 
the energy of the ground state, which would otherwise be infinitely negative. 

What we consider here is simply a free particle confined to move on a finite interval. 
Its Hamiltonian is just a free Laplacian and the only non-trivial dynamical aspect of the 
system is the behaviour of the particle at the end-points. This behaviour is described, in 
mathematical terms, by a set of boundary conditions to be imposed on the wave functions. 
The most general choice of the boundary conditions compatible with the requisite of self- 
adjointness for the Hamiltonian is parametrized by a U{2) matrix. 

Upon discretisation, the interval is replaced by points and the Laplacian by the 
classic finite difference operator already used by Laplace. In the continuum limit, one 
recovers the self-adjoint extension of the Laplacian characterized by periodic boundary 
conditions for the wave functions. One can then wonder whether it is also possible to 
recover all the other self-adjoint extensions as continuum limits of some suitably modified 
finite difference operators. Indeed, after noticing that the classical finite difference oper- 
ator is well defined only away from the end points of the interval, in order to define it 
everywhere, one is led to consider its most general symmetric extensions. The extended 
operators depend on 4 real parameters, which are the discrete analogue of the U{2) ma- 
trix found earlier in the continuum. The continuum limit is then obtained by taking 

to infinity. The corresponding four- dimensional renormalization fiow turns out to be 
exactly solvable and neverthless non trivial. We find and infinite number of fixed points 
which are related to the self-adjoint extensions of the Laplacian characterized by a set of 
scale-invariant boundary conditions. The other non-trivial trajectories, which in the limit 
approach the fixed points, correspond instead to the remaining choices of the boundary 
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conditions. 

The paper is organized as follows. Section 2 reviews some basic facts about the self 
adjiont extensions of the Laplacian on an interval. In Section 3 we introduce the finite 
difference operators which represent the discretised version of the Laplacian with arbitrary 
boundary conditions for an interval and in Section 4 we discuss the related eigenvalue 
problem. The continuum limit is discussed in Section 5. Finally, the Appendix contains a 
detailed analysis of the eigenvalue equations associated to the finite difference operators 
introduced in Section 3. 

2. Self- Adjoint extensions of the Laplacian on an interval. 



Let us consider the interval of the real line [0, 27r] and let £^([0,27r]) be the Hilbert 
space of square integrable functions over [0, 2n] with inner product given by: 

/■27r 

{(f),^) = / dx(j)*4j . (2.1) 
Jo 

Consider now the Laplacian operator: — As it stands it is only a formal differential 
operator. In order to obtain a well defined operator in £^([0, 27i]) it is necessary to specify 
the domain of functions on which it acts. As a starting domain, consider the set of 
functions which are twice different iable and vanish at the end-points together with their 
first derivatives: 

D° = {^{x) : ^(x) G C2[0,27r] ; ^(0) = V>(27r) = 9,^(0) = 9,.^(27r) = } . (2.2) 

Let us call A*^ the operator obtained by restricting the Laplacian to V^. A" is symmetric, 
ClS db simple computation shows: 

Jo 

= (AV,/) V/,V'°eP°. (2.3) 
Consider now the adjoint A°^ of A°. It is defined by the equation: 

(A°V,0°) = (V^,AV°) V/gP°. (2.4) 

The domain of A''''' is defined to be the set "D^^ C £^([0,27r]) for which equation ( ^.4[ ) is 
fulfilled, for all 0° G T>^. It is easily seen that 

P° C . (2.5) 

In fact, eq. ( p.4|) is satisfied by all twice differentiable functions ip, regardless of their 
boundary values or of the boundary values of their derivatives. This shows that A°, 
though symmetric, is not self-adjoint (s.a.). In order to make it s.a. it is necessary to 
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extend its domain from T>^ to a larger one V. If it is possible to do it and in how many 
different ways, is established using the deficiency index theorem . This theorem can be 
stated as follows. Let be a symmetric operator on some Hilbert space Ti, with domain 
T>[A^) and let A^"^ be its adjoint. Let now n_(. , n_ be the number of linearly independent 
solutions of the equation 

A^'^ijf = ±zijf , j = l,---,n± . (2.6) 

?T,_|_,n_ are called respectively positive and negative deficiency indices. Then, A^ admits 
s.a. extensions if and only if n+ = n_ = n and they can be all obtained in the following 
way. Let U = {Uij} be an arbitrary n x n unitary matrix. Consider now the domain: 



= P° span{#) + f/.,^j^^}. (2.7) 



The operators 



where 



Au : Vu^n 



A^P = + ^c,{^P)-^' - f/,,^rO, (2-8) 



V; = ^o + c,(^f)+f/,fc^[-^) 

represent, for each U, a different extension of A^ and it can be shown that they are all 
essentially s.a.. Moreover it can be proven that all possible s.a. extensions of A^ can be 
obtained in this way, by taking different unitary matrices U in equations (|2.7| ),( p78D . 

In the case of A° it is easy to check that n+ = n_ = 2. A" then admits an U{2) 
infinity of s.a. extensions, characterized by domains of the form (p.7|) . Now, it can be 
shown that, for all f/ G f/ (2) there exists a set of homogeneous linear boundary conditions 
(b.c.) which uniquely characterizes T>u: 

Vu = {t/j : afiV''(0) + af2^'(27r) + 0^3^(0) + a^^{2n) = ^ = 1, 2} . (2.9) 



In eq. ( |2.9| ), / denotes the derivative with respect to x. The set of all b.c. can be devided 
in six (in general overlapping) different classes, depending on the couple of quantities, 
chosen among ip{0),ilj{27!-),ip'{0) and ip'{27i), with respect to which ( |2.9| ) can be solved. 
We shall see that in order to get all b.c. it is in fact sufficient to consider only three cases, 
which we list hereafter. 

Case 1. aij are such that eq. (|2.9| ) can be solved with respect to 0(0) and 0(27r): 



0'(O) \ /n ^\ / 0(0) 
0'(27r); U vJ\(P{2n) 



(2.10) 



We show now that the coefficients u, v, w and ( in ( |2.10| ) have to fulfill certain conditions. 



If (|2.10| ) is to define a s.a. extension of the Laplacian it must be that ii ijj is a. function 



such that: 

(-V'",0) = (^,-0") (2.11) 
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for all (p fulfilling ( |2.10| ), 'ip also has to fulfill ( p.lO| ). Now we have: 

{-r, 0) - -0") = -r^T + r 0ir = 

= -V'*(27r)[C0(O) + V(p{2n)] + tp* {Q)[u(j){0) + W(p{27i)] + ^'*(27r)0(27r) - ^'*(O)0(O) = 
= [u^*{0) - CV^*(27r) - ^'*(O)]0(O) + [wiIj*{0) - v^*{27t) + ^'*(O)]0(27r) . (2.12) 
Since now 0(0) and 0(27r) are arbitrary, must be such that: 



^^'(0) 
^p'{2^T) 



-w* 



-c\ ( m \ 

V* J U(27r); 



(2.13) 



So, by comparing ( |2.12| ) with ( p^.lOD , we see that ip fulfills ( ^.lOj ) too if and only if: 

u = u , V = V , Q = —w 
and so the most general b.c. of this form can be written as 

, u and V real. (2-14) 



0'(O) 
0'(27r) 



u w\ f 0(0) 
-w* V J V0(27r) 



Case 2. such that ( p.9|) can be solved with respect to 0(0) and 0'(27r): 



/0'(O)\^ /n w 



0(0) 
0'(27r) 



(2.15) 



It is obvious that, if w 7^ 0, we can solve ( p^.l5[ ) with respect to 0'(27r) and the b.c. ( p^.lSj ) 
can be rewritten in the form of case 1. Then, we set ■?) = 0. A computation similar to 
that performed for case 1, shows that (|2.15| ) define s.a. extensions if and only if w = 
and u is real: 



/ 0'(O) 
V0(27r) 



u w 
w* 



0(0) 
0'(27r) 



, u real. 



Case 3. such that 



can be solved with respect to 0'(O) and 0'(27r): 



0(0) 
P(2n) 



u w\ f (p'{0) \ 
C {)jU'(2vr)j" 



(2.16) 



(2.17) 



Again a computation similar to that for case 1 shows that in order to define a s.a. extension 
the parameters u and v must be real, while ( = —w*. Then we have: 



f 0(0) 
V0(27r) 



u w 

-W* V 



0'(O) 
0'(2vr) 



u and i) real. 



(2.18) 



This set of b.c. is not already included in cases 1 and 2 if and only if 



uv + \w\ 



and u = 0. 
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We are then left with: 



(° °) (;;(»)) .... p..) 



It can be checked that the remaining possibihties lead to b.c. already included in the 
three cases we have examined. Then, the most general s.a. extension of the Laplacian on 
an interval is described by one of the above sets of b.c. 

As we have seen earlier, the space of s.a. extensions of the Laplacian on an interval 
is a four dimensional manifold, precisely U{2). Now, for the sake of simplicity and also 
because later it will allow us to draw three-dimensional pictures of the renormalization 
flow for the discretized Laplacians, we will restrict our attention to those extensions 
which are invariant under PT symmetry (P here stands for the Parity transformation 
which interchanges the end points a; = and x = 27r while T stands for Time- Reversal) . 
Now under P the wave function 6 transforms as: 



while under T we have: 
and so under PT we have: 



(P(/))(x) = 0(271 -x) (2.20) 
(r0)(x) = 0*(x) (2.21) 

4>{x) = {PT(j)){x) = 0*(27r - x) . (2.22) 



Let us check the constraints implied by the PT symmetry on the b.c. examined earlier. 
Case 1. The PT transformed wave function fulfills the b.c: 

{ 0'(27r) )~{-w* -uj[ 0(27r) ) ' ^^"^"^^ 
and so PT invariance requires u — —v. We then have b.c. of the form: 

(m]_( u w\(m\ .2 24) 

Case 2. A computation similar to that for case 1, shows that PT invariance requires 



(2.25) 



— 1. So we have: 

/0'(O) \^(u e^'\ ( 0(0) 
\<P{2n)) \e-^' 0JU'(27r) 

Case 3. PT invariance requires v = Q. Then we are left just with one set of b.c, those 
for an open line: 

0(0) = 0(27r) = 0. (2.26) 

We turn now to the eigenvalue problem for the PT invariant s.a. extensions listed 
above. The general solution to the eigenvalue equation: 

- 0" = (2.27) 
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is of the form 

(f)p = Ae'P'' + Be-'P\ (2.28) 



If p is real the function ( |2.28| ) is a plane wave solution and p represents the inverse 
wavelength, if p is imaginary ( 2.28|) represents a bound state solution. An eigenvalue 
equation for p is obtained after imposing on ( |2.28|) the b.c. ( p.24| ), (|2.25|) or (|2.26| ). We 
then get the following equations: 

Case 1. The coefficients A and B have to fulfill the following linear eqs: 

A{u + we^'^^'P - ip) + B{u + we'"'^^^ + ip) = 

A{ue'''''P + w* + ipe'^^'^) + B{ue-'''''^ + w* - ipe-'^^'P) = . (2.29) 
The consistency condition for this system is the eigenvalue eq. for p: 

{p^ + \w\^ — u'^)sin2np — 2upcos27cp — (w + w*)p = . (2.30) 

Case 2. The coefficients A and B have to fulfill now the following linear eqs.: 

A{u + ipe'^e^'^^'P - ip) + B{u - tpe'^e-''^^^ + ip) = 

A(e^2"P - e-'^) + B{e-'^^P - e''^) = (2.31) 

and the consistency condition for this system gives the following eigenvalue equation for 
p: 

2pcos2TTp + usin2Tip — 2pcos9 = . (2.32) 
Case 3. The coefficients A and B have to fulfill the linear system: 

A + B = 

^^i2np ^ 5g-i2^P = _ (2.33) 
The consistency condition now reads: 

sin27rp = . (2.34) 

Of particular interest among (|2.24|) , (|2.25| ) and ( p.26| ) are those s.a. extensions char- 
acterized by a set of scale invariant b.c. As we shall see later they will correspond to the 
fixed points of the renormalization fiow in the space of the lattice discretizations of the 
Laplacian on an interval. Among the b.c. considered in case 1, the only scale-invariant 
ones are: 

0'(O) = 0'(27r) = . (2.35) 

These are the b.c. for an open line with "free" ends. Among the b.c. under case 2, the 
scale-invariant ones are: 

0(0) = e*V(27r) 
(j)'(0) = e'V(27r) . (2.36) 



They can be interpreted as b.c. for a circle, threaded by a magnetic flux. Finally, the 
only scale-invariant b.c. of type 3 are: 

0(0) = 0(27r) = . (2.37) 

They describe an open line with completely reflecting end points. 



3. Discretized Laplacians on an interval. 

In the previous Section, we have seen that the Laplacian admits infinite self-adjoint ex- 
tensions on the interval [0, 2tt], each distinguished by a different domain Vu, characterized 
by a certain set of b.c. on the wave functions (f^.9D. 

Suppose now we discretize the interval by dividing it in N subintervals of width 2tt/N 
and centered around the points Xj = (i — |)^, i = 1, ■ ■ ■ , A^. Consider now the subspace 
Ti^ of stepwise functions x{^)) which are constant in each of these subintervals: 

97r Itt 

7i^^span{x^(a;)=x'^(x.) for {r - I)— < x < i— ^ = l,---,iV}. (3.1) 

Ti^ is an A^-dimensional subspace of -^^([0, 27r]). It is clear that, in the limit N ^ oo, Ti^ 
becomes dense in i3^([0, 27r]). We can represent each function in Ti^ as an A^-dimensional 
complex vector: {x^ i})) = ix^i^i))- 

Consider now the self-adjoint extension of the Laplacian Aper acting on periodic wave 
functions: 

0(0) = 0(27r) 0'(O) = 0'(27r) . (3.2) 

(This extension falls unders case 2, when -u = and w = 1.) As it is well known, this 
extension can be approximated with a finite difference operator acting on the space Ti^: 

(A^X^)(^) = -ZUx''i^ + 1) + X^(^ - 1) - 2x^(^)] ^ = l,■■■,N. (3.3) 



Here + 1 = 1 and = A^. is an overall normalization. Let us see how Aper 



IS 



recovered in the limit A^ oo. The eigenf unctions of A at are easily seen to be equal to: 

'N-l 



xZU) = r"^(^^ ~ ^ ) f = exp{imxj) m = 0, ±1, ■ ■ ■ , ± 



(3.4) 



where [a] is equal to a if a is integer and to a — 1/2 if a is half- integer. The corresponding 
eigenvalues are: 

A^ = -2Z^(co.^-l). (3.5) 

The eigenfunctions are seen to be stepwise approximations of plane waves. The exponent 
m can be interpreted as an inverse wavelength (in units of the interval length). Now 



7 



the discretized operators A converge to Aper in the following sense. Let us order the 
eigenstates of Aper according to their eigenvalues Ei: 



Ei < Ei. 



oo. 



(3.6) 



Let now r an arbitrary integer and consider the operators A , with N > r. Let us order 
also their eigenf unctions Xm(^) according to their eigenvalues A^: 



m 



AT. 



The claim is that the first r eigenfunctions Xi 



X 



■ ■ Xr (x) 



of A converge pointwise 
ipr (x) . Moreover it is possible to fix the 



oo, also the first r eigenvalues converge to 



to the first r eigenfunctions of Aper, i^ii^ 
constants Z^- so that, in the limit 
the Em- These statements are easily checked in this example. For instance, if we choose 
Z]y = TT", we see that: 



lim \^ 



lim 




cos- 



m 



E„ 



(3.7) 



The question that arises now is if we can recover all self-adjoint extensions of the 
Laplacian on the interval as continuum limits of finite difference operators similar to 
(|3.3| ). We have seen in Sec. 1 that the s.a. extensions of the Laplacian on an interval all 
have different domains. On the contrary, any finite difference operator we can write will 
act on the entire subspace Ti^ for every A^. It becomes then unclear how the continuum 
operators with their different domains can possibly emerge in the continuum limit. As 
we shall see in the sequel of this paper this is indeed the case. The operator p.3| ) can be 
generalized in a natural way to give a family of finite difference operators containing four 
real parameters. These parameters will be subject to a renormalization procedure having 
A^, the number of subdivisions of the interval, as control parameter and this will give rise 
to a flow in the parameter space. The corresponding trajectories will converge towards 
certain fixed points. We shall see that each of these trajectories leads in the limit to a 
different s.a. extension of the Laplacian. We shall argue that all s.a. extensions can be 
recovered by simply looking at different trajectories. 

Let us generalize the finite difference operator ( p.3|) on an interval. Notice that, on an 
interval, formula ( |3.3| ) certainly makes sense for z = 2, ■ ■ ■ , A^ — 1, but becomes ill-defined 
at the endpoints, i = 1 and i = N . If we then write A^ as a matrix acting on the vector 



[i) we get: 



N 



( * 
1 






* 


1 










(3.8) 
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In ( |3.8| ) the first and the last raws of the matrix representing are missing. Since there 
is no natural way of defining at the end points, we will simply take for A^ the most 
general hermitian completion of its corresponding matrix and write: 



A 



N 



y2 



1 a 


1 





1 


-2 


1 





1 


-2 












1 



c 





V c* 



1 





-2 
1 



(3.9) 



Here a and b are arbitrary real parameters, while c ia an arbitrary complex number. 
Notice that the matrix elements containing the parametrs only relate the wavefunction 
at the endpoints and seem therefore related to boundary conditions. We also remark 
that a, b, c constitute a set of four real paramters, which is also the dimension of the set 
of U{2) matrices that parametrize the s.a. extensions of the Laplacian. 



4. The eigenvalue problem for A 



N 



In this Section we study the eigenvalue equation for the matrix A (for simplicity of 
notation we will omit from now on the superscript from the symbols for the finite 
difference operators A^ or the wavefunctions x^{^))'- 



/ a 


1 


■■• 




c 


\ 


( Xa(1) \ 


1 


-2 


1 




... 









1 


-2 1 





... 











... 


1 


-2 1 






\c* 










1 b 


J 


vxA(iv); 



/ Xa(1) \ 



(4.1) 



VXA(iV)/ 



The eigenfunctions of A fall into 5 classes, depending on the eigenvalue A. They are listed 
hereafter: 

I) If A < 0, we have: 

Xx{n) = Ae^"" + Ee"^", (4.2) 



where A = —2{coshk — 1) , A; > 0. 
II) If A = 0, we have: 



Xo(^) = A + Bn 



(4.3) 
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Ill) If < A < 4, we have: 



Xx{n) = Ae'"'' + Ee-*'^", (4.4) 

where A = —2{cosk — 1) , < k < n. 

IV) If A = 4, we have: 

X4(n) = (-l)"(A + 5n) . (4.5) 

V) If 4 < A < oo, we have: 

Xxin) = AC + Br'', (4.6) 
where A = -(^ + - 2) , ^ < -1. 

We prove the statement for case V, the proof being similar in the other cases. Let xx 
be an eigenfunction with eigenvalue A > 4. Let then ^ be the unique real number smaller 
than -1 such that A = — (^ + — 2) and define A and B to be such that: 

Xxil) = A^ + Br' 

Xxi2)=Ae + Br'. (4.7) 
Notice that eq.ns ( [4.7] ) have a unique solution because 

By construction, xx is of the form ( [4.6|) for n = 1,2. We now prove by induction that it 
must be so for all n. So, suppose the statement is true for n <n, with n> 2. Then, since 
Xx is an eigenfunction of A^, it must be that: 

-[Xx{n + 1) + Xx{n - 1) - 2xxin)] = Xxxin) 



or 



Xx{n + 1) = (e + - 2){AC + BCl - {AC-' + i?r^^"'^) + 2{AC + BCl = 

= AC^^ + 5r^"+^^ , (4.8) 

which proves the statement. 

The physical meaning of the solutions of type I, II and III is clear: they represent 
stepwise approximations of bound states, zero modes and plane waves respectively. As 
for the solutions of type IV and V we shall see in the next section that they "disappear" 
from the spectrum in the continuum limit and are unphysical. Notice also that the 
eigenf unctions of type I, III and V can all be written as: 

XA(n) = + A 7^ 0,4 (4.9) 

where 

A = -{z + z'^-2), 
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z eC =]- oo,-l[\J {z = e'\ < ^ < tt} y ]1, oo[. 

Let us derive now the eigenvalue equation for the states of type I, III and V. Notice 
first the the wavefunctions ( |4.9| ) automatically satisfy the eigenvalue equation for n = 
2, ■ ■ ■ , N — 1. Then we just have to consider the endpoints n = 1 and n = N . We find : 



Ajv[(a + 2)z + cz^ - 1] + 5jv[(a + 2)z"^ + cz"'"' - 1] = 



AM[{h + 2)z^ + c*z- z^ 



z 



-(iV+l)i 



0. 



(4.10) 



The condition of compatibility for these eqs. leads to the eigenvalue equation for z: 



\a[z 



-{iV-l)^ 



+ P{z 



N 



-N 



z 



[Z 



N+1 



-{Af+l)^ 







(4.11) 



where we have defined: 



a 



|c|2- (a + 2)(6 + 2) 
/3 = a + 6 + 4 
7 = 2(Rec) . 



(4.12) 



The real parameters and 7 introduced here are independent on each other and are 
only subject to the condition: 

4a + /^^ - 72 > . (4.13) 

For the sake of simplicity we will restrict our analysis, as we did in the continous case, 
only to the set of matrices A which are left invariant by a PT transformation. Now, under 
PT the vector x(^) transforms as: 



{PTx){n) = x\N -n + l) 
and then in order for the matrix A to be PT invariant it must be that a 



A 



(4.14) 



h: 



1 a 


1 









c 




1 


-2 


1 







... 







1 


-2 


1 





... 














1 


-2 1 
















1 a 


) 



(4.15) 



A PT invariant A is thus identified by three real constants and we can think of it as a 
point in R^. In the next Section we will turn to the study of the continuum limit for the 
finite difference operators ([4.15|) and show that in this limit we recover all PT invariant 
s.a. extensions of the Laplacian, eqs. (|2.24| ), (|2.25| ) and (p.26|) . 
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5. The continuum limit for A . 



In this Section we study the continuum hmit of the finite difference operators ( ^4.15| ). 
In the Appendix it is shown that at least N — 2 eigenstates of A are of the plane- wave 
type (eq. (^4.4D). Let us rewrite them: 



Xfc(n) = v4e*'='^ + 5e-^'=" = Ae^^" + 5e-*^'^ = XpH, < A; < vr . (5.1) 

The number p = ^ represents the inverse wavelength of the state in units of the interval 
length. In the Appendix it is shown that the eigenvalues k are distributed essentially 
uniformly in the interval ]0, 7r[, which implies that the first values of p are of the order of 
1, 2, ■ ■ -. The lowest eigenstates then have wavelengths of the order of the length of the 
interval. 

Now, our prescription for taking the continuum limit is the following. Assume we start 
with a certain set of values for a, c and N , (ao, Cq; Nq), or equivalently with a set of values 
for a, (3, 7 and A^, (ao, Po, 7o; A'o)- We then consider the eigenvalue equation determining 
the wavelengths of the "plane wave" states ( [4.4|) , namely eq. (|4.11| ) with z = e*^'^^/^: 



aosin2Tcp^ — —j + i3osin2TTp + 7o,sm27rp— = sm27rp^l + —j . (5.2) 

Let now: 

Pi<pI<---<P- ^<N, (5.3) 

be its solutions (here we just list the eigenvalues and do not consider the existance of 
eventual degeneracies, which, as it is shown in the Appendix is an exceptional case). Our 
prescription for taking the continuum limit is then to vary a, (3 and 7 with A^ so that the 
corresponding eigenvalue equation for the plane wave states: 

a{N)sin2'np{l - + l3{N)sin2-np + 'y{N)sin2TTp^ = sin2TTp(^l + (5.4) 



shares with ( p.2| ) as many wavelengths as possible, starting from the longest one. We shall 
see that this prescription fixes completely the dependence of a, f3 and 7 on A^. In general, 
it will be possible to hold constant only the first three inverse wavelengts, p^, P2, p^, but 
there will be exceptional cases where we will have to consider the first four of them. 

Having determined the evolution of the parameters a,/3,7 with A^, we have to assign 
now the renormalization condition for the overall constant Z^v which appears in front of 
A^, eq. ( [4.151 ). This constant can be fixed by demanding that the eigenvalue Ei = ZffXi 
of the first plane-wave state (the one with inverse wavelength p^) takes some fixed value, 
independent on A^. Now, it can be verified, by studying the eigenvalue equation, that 
Pi is always of the order of unity, which implies for ki a value of the order of 1/A^. 
Consequently, we have: 

E, = ZjjXi = -2ZUcosk, - 1) ^ {Z^kif = (zn"^] ■ (5.5) 
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We see from this formula that Z\ diverges in the continuum hmit hke A^. Notice that 
the scahng behaviour of imphes that the eigenvalues of the states of type IV and V, 
which are always larger than or equal to 2^^4, diverge in the limit N oo and disappear 
from the spectrum. This means that the states of type IV and V do not have physical 
meaning. 



Our renormalization procedure generates a flow in the space of PT invariant matrices 
A^, the control parameter being A^. We will then be able to prove that, depending on 
the choice of the initial values Oq, (^o, 70) the discretized Laplacians ( [4.15 ) converge in the 



continuum limit to any of the PT invariant s.a. extensions discussed in Sec. 1. This will 
be achieved in two steps. First we will prove that the eigenvalue equations (|5.4|) converge 
in the limit ^ oo to those found for the PT invariant s.a. of the Laplacian, eqs. 
(|2.3CI|) , ( |2.32|) and ( |2.34| ). The wavelength spectra of the finite difference operators 
then converge to those of the s.a. extensions of the Laplacian. Then we prove that also 
the eigenfunctions converge, in the sense discussed earlier, to those of the continous 



PT invariant Laplacians. This will be done by showing that the linear systems ( [4.10|) 



which determine An and Bn converge to the corresponding systems for the coefficients 
A and B of the continuum eigenfunctions. 

We summarize here the qualitative features of the flow and anticipate what boundary 
conditions for the Laplace operator one gets in the continuum limit: 

• We can describe the space of PT-invariant matrices A^ by means of the three 
parameters a, b , c or, alternatively, by the real numbers a, f3, 7 defined in ( [4.12 ) 
and subject to condition ( [1.13| ). Hence the parameter space is: 

V = {{a, (3, 7) G R3 : 4a + /?2 _ ^2 > Q| _ 

• The renormalization group has an infinite number of fixed points, belonging to the 
sets: 

Li = (a = -l,|/3|>2,7 = 0) 
L2 = (a = l,/3 = 0,|7| <2) . 

• If the initial data ao, Po, 70 are chosen to belong to such lines, in the continuum limit, 
one recovers the s.a. extensions of the Laplacian on an interval which correspond 
to scale-invariant boundary conditions. In particular the point P = (—1,2,0) G Li 



corresponds to the s.a. extension (|2.35[ ), the other points on Li to the s.a. extensions 



( lO^ ) and the points on L2 to ( [X^ ) 



For all initial data not belonging to a certain surface S{ao, Pq, 70, Nq) = OmV, the 
parameters a, (3, 7 follow trajectories that all approach the point P, for N 00 
(Figure 1). The continuum limit Laplace operators corresponding to such trajec- 
tories are the ones characterized by the non-scale invariant boundary conditions 
falling under Case 1. 



13 



• If the initial data are chosen to belong to the surface S, the trajectories approach 
a point on the line L2 (Figure 2). The corresponding continuum limit Laplace 
operators are charcterized by the non-scale invariant boundary conditions falling 
under Case 2. 

Let us show now that our renormalization prescription is well posed and prove the 
above statements. 

First, consider the case when ao,/3o,7o belong to Li. It is easy to verify that the 
eigenvalue equation ( [5 .21 ) reduces to 

N , , 

sin2'Kp = < p < — . (5.6) 

The eigenvalues are the integer and half- integer numbers m less than A^/2. Consequently 
the only effect of increasing is that of adding shorter wavelengths to the spectrum, 
whithout changing the preexisting ones: according to our prescription these points then 
represent fixed points of the flow. Let us look now at the behaviour of the eigenfunctions 
Xm- If /^o = 2, we can easily verify, by expanding z = exp{i2Trm/N) in powers of 1/N 
and then keeping only the leading terms in eqs. (|4.10|) , that these equations reduce 
to those determining the coefficients A and B of the continuum eigenfunctions for the 
Laplacian with b.c. given by eq. (|2.35|) . We can then conclude that the operators 
associated to the fixed point P = (—1,2,0) converge to the extension ( |2.35| ). If 
/5o 7^ 2, a similar computation shows that in the continuum limit the eigenfunctions of 
the discretized operators converge to those for the s.a. extension of the Laplacian with 
the b.c. given by eqs. (|2.37|) . 



Suppose now that the initial point belongs to L2. Eq. (|5.2|) reduces to: 

N , , 

2cos27rp = 7o; < p < -— . (5.7) 

Again the only effect of increasing is just that of adding shorter wavelengths to the 
spectrum and then all points of L2 are fixed points. As to the eigenfunctions Xj^; by 



looking at the limit of eqs. ( [4.10|) , we can check that they converge to those for the s.a. 
extension of the Laplacian having the b.c. (|2.36|) , with IcosQ = 70. 

We consider now the general case. Let pi, p2 and ps denote the first three eigenvalues 



of eq. (|5.2|) and suppose that the initial values of the parameters ao, /?o, 7o; -^0 are such 



that the entire function of the complex plane: 

h{x;p^) = det || sm27rp°(l — x), sm27rp°, sin27cp^x \\ j = 1,2,3 (5.8) 

does not vanish identically in x [h is identically zero if and only if at least two of the 
three eigenvalues p?,P2?P3 are both integer or half-integer; this exceptional case will be 
examined later]. Since now h{x;p^) is an entire function of x, its zeros will be isolated 
points in the complex plane and this implies that there will be at most a finite number 
of them in the disc |x| < 1/A^o- Now, according to our renormalization prescription, the 
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running constants a{N), j3{N) and 7(A^) will be chosen to be the solutions of the system 
of three linear equations: 

a{N)sin27rp^j ^1 - + /5(A^)sm27rp° + 7(A^)sm27rp°^ = sm27rp° + . (5.9) 

This system admits a unique solution for all N > Nq [except possibly for a finite number 
of values of such that is a zero of h{x). When this is the case we simply neglect 
these values of A^]. This proves that our scheme is well posed. 

We will now analize the limiting behaviour of the solutions a{N), j3{N) and 'j{N). For 
this purpose, we define two more entire functions of x: 

f{x;p^) = det II sm27rp°(l — x), sm27rp°(l + x), sin2iTp'^x \\ (5.10) 

and 

g{x]p^) = det II sm27rp°(l — a;), sm27rp°, sm27rp°(l + x) ||; (5-11) 

notice that /(x) is an odd function of x, while g{x) is even and that /, g and h all vanish 
at a; = 0. The solution of the system ( p.9|) can now be written, in terms of /, g and h as: 

/i(-x;p°) 

«(-;P.) = -l^^> (5-12) 

n[x;pj) 
h{x;pj) 

a{x;p'j),P{x]p'j) and 7(0;; p^), being ratios of entire functions, are thus analytic functions 
of the complex plane having at most poles of finite order. We are now interested in the 
limits of these functions when x —>■ 0. Let us start from a{x]p^). Let n{p^) be the order 
of the first non vanishing derivative of h{x\p^) for x = (fi represents also the order of 
the zero of h{x) at the origin; a simple computation shows that n(p°) is always larger 
than 1). Use of the I'Hopital theorem then implies: 



The limiting value of a is then always equal to either 1 or -1. In order to study the limits 
of /? and 7 let us write the power expansions of f and g around x = 0: 

00 

h{x;p'^)= K{p])x^ (5.16) 

n=n(p°) 
00 

/(x;p°) = E/2..(P,V" (5.17) 

n=l 
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0^ ™2n+l 



jX 



(5.18) 



n=0 



where hn, fn, Qn denote the derivative of order n with respect to x of h, f , g respecively. 

In order to recover, in the continuum hmit, all s.a. extensions of the Laplacian, it 
seems sufficient to consider only two cases n{p^) = 2, 3. 

Case a. The initial data are such that n{p^) = 2. This is the most general case, h then 
has a zero of the second order at the origin. According to (|5.15| ) we then have: 



lim a(x: n' 



-1 . 



(5.19) 



It is easy to check now that: 



and then use of the I'Hopital formula in eq. ( |5.13| ) implies: 

/2(P°) 



2. 



(5.20) 



As for 7(x) it can be checked that (72 (Pj) = and then we find: 

x^o J h2{Pj) 

By putting together eqs. ( ^.191) , ( |5.20|) and (|5.21| ), we see that: 

if n{p'^) = 2 lim(«(x),/3(x),7(x)) = (-1,2,0) = P. 



(5.21) 



(5.22) 



We saw earlier that the limiting point P is a fixed point of the fiow and that the discretized 
operator associated with it converges in the continuum limit to the s.a. extension of the 
Laplacian characterized by the b.c. ( p.35|) . 

We show now that the eigenvalue equation ( |5.4| ) converges, for iV — >■ 00, to an equation 
of the form (|2.30|) . For this purpose, let us expand the r.h.s. of eqs. ( ^.12|) , ( p.l3| ) and 
( |5.14| ) in power series of x = 1/A^: 



a(x;p°) = -1 + i/,(p°)l + H^ip'^) (^1^ + 
/3(x;p;) = 2 + F,(p;)l + F,(p;)(iy + .. 



(5.23) 

(5.24) 



7(x;p;) = Gi(p;)l + G2(p;)(^) +••• . (5.25) 
If we substitue these expansions in ( p. 4] ) and keep only the leading terms in 1/A^ we get: 

( - 2 + H,{p^)^ + i^2(p?)^)co.^ + 2 + F,(p;)l + F,(p;)i, 



sin2'Kp 



16 



(5.26) 



N 

A simple computation shows that Hiijf-) + Fi{p^) = and then eq. (|5.26|) reduces to: 

[{2npy + if2(p?) + F2(p°)]szn27rj9 - 27ip[Hi{p''j)cos27ip - ^(pj)] = (5.27) 
This eq. has the same form as (|2.30|) , if we identify: 



1 

\wr — u 



(27r)2^ -v^^ -^n' 
1 

2^' 



2m = — ifi(p5 



w + w* = -^Gi{p'l) . (5.28) 

Once these identifications are made, it is possible to check, by susbstituting the expansions 
(|5.23| ), ( |5.24| ) and ( [5.25| ) into eqs. ( [4.10| ), that we recover, to the leading order in 
eqs. ( p.29| ). This shows then that the continuum limit of our discrete operators coincides 
with the PT invariant s.a. extensions of the Laplacian contemplated under case 1. 

By varying ao,/?o,7o or equivalently p° we can get essentially all s.a. extensions of 
type 1. Consider infact eq. ( [^.301 ). Our prescription for the continuum limit is such that 
2;g]2iain eigenvalues also in that limit and then they are roots of eq. ( p.30|) , with u and 
w given by eqs. ( p.28| ). Now, conversely, it can be checked that, if we start from eq. ( p^.30 ) 
and try determine the coefficients u and w such that p° are solutions, the answer is in 
general unique and is given again by eqs. (|5.28| ). So, assuming that by varying ac/^cTo 
we can get all the triples p° that can be obtained from eq. ( |2.30| ) by looking at its first 
three solutions, we can conclude that all of s.a. extensions of type 1 can be recovered. 
There are a only few exceptional cases that are left: they correspond to the case when two 
or more among the are integers or half-integers. These last extensions will be found 
by analizing the case of ao,/3o,7o such that h{x;p^) = 0. 

Case b. Suppose that the imitial data (oo, /3o, 7o; -/Vq) are such that n{p'j) = 3. This is 
the same as saying that: 

/i2(p°) = det II -2np°jCos2np°j, sin2np°j, 2Txp] \\= (5.29) 

and 

h^{p]) = -det II -{2Tip']fsin2'Kp], sin2Trp], 27rp° ||^ . (5.30) 

These eqs. specify a surface S = ^(ao, /5o, 7o; iVo) = in the parameter space. Now eq. 
(|5.29|) implies then that there exist three real numbers /i, u and p such that: 

/x27rp°cos27rp° + z/sm27rp° + p27rp° = . (5.31) 

Moreover eq. ( |5.30| ) implies that fi ^ 0. h{x) now has a zero of the third order at x = 
and eq. ( ^.15|) implies that: 

lima(x;p°) = 1 . (5.32) 
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A simple computation also shows that: 



lim 3(x: p. 



lim 7(x; p' 



fx 



(5.33) 
(5.34) 



Condition ( [4.13| ) then implies that |p/yu| < 1. 

This computation shows that the trajectories starting from a point of 5* will approach 
the point (a,/3, 7) = (1,0, — 2p//i) belonging to the set L2. We showed earlier that the 
points of this set represent fixed points and that in the continuum limit they correspond 
to the s.a. extensions of type 2, characterized by the b.c. ( [^.3(j| ). If we now expand the 
r.h.s. of eqs. ( p.l2| ), ( |5.13|) , ( ^.14|) around x = we find: 



/?(a:;p,") = Fi(p^")- + F2(p; 



N 



N ^ 
2 



+ ■ 



7(x;p° 



(5.35) 

(5.36) 
(5.37) 



Notice that these expansions do not coincide with those conputed earlier, under case 1. 
If we substitute these expansions in the eigenvalue eq. (|5.4| ) and keep only the leading 
order terms in 1/iV we get the limiting eigenvalue eq.: 



cos2ttp + — 







(5.38) 



As a consequence of eq. ( |5.31| ), we have the following relations: 

2u 



Giip') 



3/x 

V 

3/i 
2pv 



and then eq. ( |5.38| ) can be rewritten as: 



— sin2'Kp + Airp 
/i 



cos2'n-p + — 



which is of the same form as eq. (|2.32 ), if we identify: 



u 



27r/x 
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(5.39) 



(5.40) 



(5.41) 



cose = -- . (5.42) 

Using these identifications, we will now prove that, in the continuum limit, the linear 
system ( [4.10| ) is equivalent to its analogue, eqs. ( p.31| ), for the PT invariant extensions 



contemplated under case 2. Notice first that, to the zeroth order in l/N, eqs. (|4.1CI|) are 
both equivalent to the following eq.: 



^(^^i0^i27.p _ ^) ^ 5(e*^e-*2"P - 1) = (5.43) 



which is clearly equivalent to the second of eqs. ( |2.31| ). If we now subtract the second of 
eqs. ( [4.10| ) from the first and expand the difference to the first order in and use eqs. 
( ^ ), we find: 

+ ipe-'-^e'^^P - ip \ +B\^ - ipe-'-^e-'^''^ + zp] = (5.44) 
\27r/i J \2nfi J 



which, in view of eqs. ( |5.41| ) and (|5.42| ), is equivalent to the first of eqs. ( |2.31| ) 



We now turn to the exceptional case, when the initial data {ao, Po,'yo; Nq) are such 
that h{x;p^) = 0. This case occurs if and only if at least two of the eigenvalues p^ are 
both integer or half-integer and this happens if and only if a + 1 + 7 = 0. In this case, 
Pi and P3 are actually half-integer, while P2 is neither integer nor half-integer. The plane 
of equation a + l + 7 = 0is then invariant under the renormalization flow , but in order 
to completely fix the evolution of the running constants <y{N), (3{N) and 7(iV) we have 
to consider also the fourth eigenvalue p^, which again is neither integer nor half- integer. 
If we use the relation a + l + 7 = 0to eliminate 7, the evolution of a and /3 is then 
determined by the equations: 

a{N)[sin2TTp^ {1 — x) — sin2np^x] + P{N)sin2np^ = sm2Txp^{l + x) + sin2np^x I = 2,4: 

(5.45) 

(in what follows we will not show again the range of values for the index /; it will be 
understood that / = 2, 4). We can apply to these eqs. the same methods as those used for 
eqs. ( ^.4|) , and the conclusions are similar. Therefore we will be brief here. We introduce 



the entire functions: 

s{x]p^) = det II sin2iTp^(l — x) — sin27rp^x, sin2TTp^ \\ (5.46) 

t(x;p°) = det II sm27rp°(l — x) — sin2Txp^iX, sin2'np^{l + x) + sin2T[p^x \\ . (5.47) 

Notice that both these functions vanish for x = and that t(x; p^) is an odd function of 
X. Using these functions we can write the solution to eqs. (|5.45| ) as: 
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Let now r(p°) be the order of the first non-vanishing derivative of s{x) for x = 0. Use of 
I'Hopital theorem implies: 

Vima{x;p'}) = {-iy^p"l (5.50) 
As before, it is sufficient to consider two cases. 

Case a: r{p^) = 1. A simple computation analogous to that done earlier shows that: 

Iim(a(x;p0),/3(a;;p0),7(x;p°)) = (-1,2,0) . (5.51) 

Moreover, the eigenvalue equation, eq. ( [4.11| ) and the linear system determining A and 
B , eqs. ([4.10|) , reduce in the continuum limit to the corresponding equations for the 
Laplacian with b.c. of type 1. 
Case b: r{pf) = 2. This case occurs if and only if: 



ds{x]p^i) 



= —det II 2np^{cos2np^ + 1), sin27rp^ 

x=0 



dx 

— 4:Cosnp2COSTcp^det \\ 2'rrp^cosTcp^, sirnvp^ 11= . (5.52) 



This conditions is fulfilled if and only if there exist real numbers a and r such that: 

27iap^sin7rp^ + rcosnp^ = . (5.53) 
Computations similar to the previous ones show that: 

lim>(x;p°),/?(x;p°),7(x;p?)) = (1,0,-2) . (5.54) 



The limits of the eigenvalue eq. (|4.11|) and of the linear system ( [4.10|) coincide now with 
those for the Laplacian with b.c. of type 2. 
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APPENDIX 



A. Study of the eigenvalue equation for 

In this Appendix we analize the eigenvalue equation for the finite difference operators 
A^. We start by showing that at least N — 2 eigenf unctions are of type III, (eq. ( [4.4| )), 
namely of the plane- wave type. Consider the eigenvalue eq. ( [4.11|) , for z = e^^ with 
< < tt: 

asink{N - 1) + (3sinkN + '^sink - sink{N + 1) = . (A.l) 

< A; < TT . 

It is convenient to expand the sines in eq. ( |A.1|) and rewrite it as: 



[{a - l)cosk + p]sinkN = [{a + l)coskN - 'j]sink . (A.2) 

There are now several cases. Here we will examine only some of them, the analysis being 
similar in the other cases. First, there are some exceptional cases. 

Case a: a = —1 and 7 = 0. Condition ( [4.13| ) then implies \f3\ > 2. These points all 
belong to the line Li. In this exceptional case eq. ( [A. 1|) simplifies to 

sinNk = sin2TTp = . (A. 3) 

The eigenvalues p are the integers and the half-integers less then N/ 2 and there are N — 1 
of them. It can be checked that N/2 is also an eigenvalue. 

Case b: a = 1 and f3 = 0. Condition ([4.13|) implies I7I < 2. These points all belong 
to the line L2. In this other exceptional case eq. ( |A.1|) reduces to: 

2coskN = cos27cp = 7 . (A. 4) 

Case c: a + l±7 = 0; a — 1 + (3 ^ and (3 and 7 are not such that we are in 
cases (a) and ( b). Now kN = {2m — l)7r,m = 1, ■ ■ ■ [A/2] are solutions to ( |A.2| ). The 
remaining solutions are the roots of the equation: 

(a-l)cosk + /3 , ^, , kN ^. 
^ ' - (a + \)cotg . (A.5) 



sink 2 

This equation always has exactly one solution in each interval ;^(2m — 1) < /c < j^{2m + 
1), m = 1, ■ ■ ■ , [A/2] for A odd and m = 1, ■ ■ ■ , A/2 — 1 for A even. These solutions are 
neither integer nor half-integer. All in all we have then A — 1 plane wave solutions. 

Case d: a + l±77^0 (and the values of a, (3 and 7 are such that we are not in case 
(b). This is the most general case. Now p cannot be integer or half-integer and then eq. 
( [A.l| ) is equivalent to: 

(a — l)cosk + 13 (a + l)coskN — 7 (A 6) 



sink sinkN 
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There are now several subcases, depending on the relative signs of a — 1 ±/5 and a + 1 ±7. 
We will discuss here only two subcases, the analysis being similar in the others. So suppose 
that: 

;9>|a-l| |7|<a + l. (A.7) 

Condition (|4.13|) is then always satisfied. The l.h.s. of eq. (|A.2|) is then always positive 
for < k < 71 and diverges both at k = and k = tt. The r.h.s. is a periodic function 
with period 2t[/N . Its qualitative behaviour is the same as that of the function cotgkN 
and it is easy to verify that there is exactly one solution in each of the intervals < 
k < ^{n + 1), n = 1, ■ ■ ■ , N — 1. There are then at least — 1 plane wave solutions. 

It is easily verified that there is also one solution of type I, namely a bound state. Since 
we have managed to find A^ eigenvalues, it also follows that they are all non degenerate. 

Consider now the case: 

(3> \a-l\ 7 > |a + 1| , (A.8) 

where a, (3 and 7 have to fulfill condition ( ^.13| ). The l.h.s. has the same behaviour 
as in the previous subcase, but the r.h.s. has the qualitative behaviour of the function 
—'-f/sinkN. It can be checked that the conditions ( 4.13|) ensures that the minimum value 



of the l.h.s. is larger than the minimum value of the r.h.s.. It is then easy to verify that 
there are exactly two distinct eigenvalues in each of the subintervals jf{2m — 1) < k < 
jj'2m, m = 1, ■ ■ ■ [A^/2] — 1. Moreover, if A^ is even there is one more solution in the 
interval ;^(A^ — 1) < A; < For all A^ there are then A^ — 1 solutions. 

Again, it can be checked that there exists one solution of type I. As before, these 
eigenvalues are then all non-degenerate. 



We show now that eq. (|A.2|) admits degenerate solutions only if a,/?, 7 belong to a 
a certain curve. From eqs. ( |4.1CI| ),with a = b and z = e^^, we see that A; is a degenerate 



eigenvalue if and only if the coefficients oi An and Bn all vanish simultaneously. This 
condition is equivalent to: 

(a + 2)e''' + ce'''^ = 1 
(a + 2)6-"" + ce-'"^ = 1 . (A.9) 
This system is compatible if and only if sin{N — l)k 7^ 0. Then its solution is: 

/? sinNk 



a + 2 



2 sink{N 
sink 



(A.IO) 

sink{N - 1) ^ ' 

When k is varied between and vr in eqs. (|A.10|) 7 describe a curve: this is then 
the set of values of the parameters for which some of the eigenvalues of eq. ( |A.1|) can be 
degenerate. 
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Fig.l. Some examples of trajectories approaching the point P. 



